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We study the pionic decay of a possible dibaryon d! — > N + N + ir in the 
microscopic quark shell model. The initial d! dibaryon wave function (J P =CP, 
T=0) consists of one 1%uj six-quark shell-model s 5 p[51]x configuration. The 
most important final six-quark configurations s 6 [6]x, s 4 £> 2 [42]x and (s 4 p 2 — 
s 5 2s)[6]x are properly projected onto the NN channel. The final state NN 
interaction is investigated by means of two phase-equivalent - but off-shell 
different - potential models. We demonstrate that the decay width depends 
strongly on the short-range behavior of the NN wave function. In addition, 
the width T^' is very sensitive to the mass and size of the d' dibaryon. For 
dibaryon masses slightly above the experimentally suggested value M<j' =2.065 
GeV, we obtain a pionic decay width of IV ~ 0.18-0.32 MeV close to the 
experimental value ~ 0.5 MeV. 
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I. INTRODUCTION 



During the last decade much attention has been devoted to theoretical and experimental 
investigations of the pionic double charge exchange (7rDCX) process on nuclei. Because this 
reaction tt + + (A, Z) — > (A, Z + 2) + tt~ involves (at least) two nucleons in the nucleus, the 
7rDCX cross section depends sensitively on short-range NN-correlations in nuclei. There- 
fore, it provides a good testing ground for the nucleon-nucleon interaction at short range. 
Experiments on different nuclear targets have unambiguously confirmed the existence of a 
narrow resonance- like structure in the 7rDCX cross-section at small incident pion energies 
T n ~ 50 MeV [|J. The position of this peak turned out to be largely independent of the 
studied nucleus. The height and width of this peak could not be explained by standard 
calculations based on the two-step process || (n + n + 7r + — > n + p + 7r° — > p + p + 
So far, these data could only be explained with the assumption of a non-nucleonic reaction 
mechanism proceeding via an intermediate dibaryon resonance, henceforth called d! . 
The quantum numbers of the d! dibaryon candidate were determined as J p =0 _ , T=0, and 
its free mass and hadronic decay width were suggested to be Md'=2.065 GeV and IV — 
0.5 MeVQ. More than a decade ago Mulders et al. || predicted a dibaryon resonance with 
quantum numbers J p =0~, T=0 and a mass M ~ 2100 MeV within the MIT bag model. 
Recently, this dibaryon candidate has been investigated in a series of works |§-|3| within the 
Tubingen chiral constituent quark model. These works emphasize the crucial role of the 
confinement mechanism for the existence of the d! . 

The quantum numbers J p =0~, T=0 of the d! resonance prevent the decay into two 
nucleons and the only allowed hadronic decay channel of the d! is the three-body decay 
into a 7rNN system with S-waves in each particle-pair 0,0] • Because the dl mass is 
only m 50 MeV above the ttNN threshold, the d' decay width IV should be anomalously 
small owing to a very small phase volume of three-particle final states. We recall that the 
currently available experimental evidence of dibaryon excitations in nuclei is very limited 



This value is uncertain by a factor of two Q . 
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0. This is due to very large N-N decay widths of most dibaryon resonances, which renders 
them undetectable on the background of other hadronic processes at intermediate energy. At 
present, the experimental evidence for narrow dibaryons is reduced to a single candidate, the 
g?'(2065). In contrast to the deuteron, which consists of two on the average widely separated 
nucleons, there are indications ||[7]], that the d! is a rather pure, compound six-quark system. 
Therefore, the dynamics of its hadronic decay into the 7rNN system should be sensitive to 
the overlap region of the two outgoing nucleons; a situation that is ideal for understanding 
the role of quark degrees of freedom in the short-range nucleon-nucleon interaction (see e.g. 



Ref. |T(| and references therein). 



Starting from this point (for alternative approaches see in Refs. [p|, pTT] , [T2[] ) we consider 
the d' decay as a (quark) shell-model transition from one six-quark configuration to another 
one by emitting a pion. The quark line diagram of the decay is sketched in Figure [I]. The 
calculation of the transition matrix elements d! — > N + iV + tt is similar to the calculation of 
A-isobar-decay matrix elements A — ► N + tt (spin and isospin flip of a quark). In the case of 
the d' decay only the initial dibaryon state is a definite six-quark configuration (the lowest 
shell-model configuration with quantum numbers J p =0~, T=0), whereas the final state 
consists of a continuum of NN-states which have to be projected onto a basis of six-quark 
configurations with quantum numbers J p =0 + , T=l of the NN 1 Sq wave. The main difficulty 
in comparing the calculated width with experimental data is its sharp dependence on the 
energy gap between M# and the 7rNN threshold. A reliable result on can be obtained 
only if the exact mass in vacuum is known (e.g. from electroexcitation of the d' on the 
deuteron at large momentum transfers [0). At present, we have only indirect data in the 
nuclear medium [[[]]. Due to the absence of vacuum data, we investigate the problem of the 
d! decay width starting from theoretical quark-model results ||,[7| for and the hadronic 
d' size parameter 6 6 . 

Our first calculation for IV was published in Ref. ||. The aim of the present work is 
to improve mainly on three important effects which were neglected in Ref. ||: a) antisym- 
metrization of the final NN-state on the quark level taking into account the effect of quark 
exchange between the two nucleons at short range, b) insertion of a complete basis of fi- 



nal six-quark states including besides the non-excited s 6 shell-model state all Pauli-allowed 
excited configurations s A p 2 and s 5 2s, which have a non-vanishing overlap with the final NN- 
state and can be populated via the emission of the pion from the initial d! dibaryon, and c) 
inclusion of the final state interaction (f.s.i.) for the two-nucleon system. 



II. DECAY DYNAMICS IN TERMS OF QUARK DEGREES OF FREEDOM 

A. Initial state 

As in Ref. || we consider only the simplest six-quark configuration s 5 p[51]x in the initial 
state (the energetically lowest J p =0~, T=0 translationally invariant shell model (TISM) 
state which satisfies the Pauli exclusion principle). It has been shown in that the d' 
wave function may be considered as a compound six-quark state, for which a single shell 
model vector provides an adequate description. This state vector is defined by 

|<0 = \s 5 p(b 6 ) [51] x , [2 3 ] C [3 2 ] T ([2 2 1 2 ] CT )[42] 5 : [21% TS} LST =110, J p = 0") . (1) 

The characteristic oscillator parameter in the six-quark wave function b§ may for example be 
determined from the minimisation of the d' mass for a given microscopic quark-quark Hamil- 
tonian The Young schemes [/d], D=X,C,S,T in orbital, color, spin and isospin space, 

as well as for the coupled spaces CT, CTS, are necessary for the unambigious classifica- 
tion of shell-model basis vectors in terms of irreducible representations (IR) of the following 
reduction chain for unitary groups: 

SU(2A) XC st D SU(2) X x SU(12) CST D SU(2) X x SU(Q) CT x SU(2) S D 

SU(2) X x SU(3) C x SU(2) T x SU(2) s (2) 



The fractional parentage coefficient (f.p.c.) technique |14-[17| based on scalar factors (SF) 



of Clebsch-Gordan coefficients of the above group [|16| -|T9|| , sketched in the following section, 
is used for the calculation of matrix elements and overlap integrals. 



B. Transition operator 



The pionic decay width of the d! is calculated, as in Ref. ||, assuming a direct coupling 
of constituent quarks with the isotriplet of pion fields <fi through the operator 

^^g^-Kr^) -"-^^' . (3) 
Here, Tj, crj and Tj are coordinate, spin and isospin of the j-th quark, k is the pion momentum 



in the c.m.s. of the d', and E n = ym 2 + k 2 . f wq is the irqq coupling constant. Its value is 

f 2 

connected with the ttNN coupling f^^ (we use = 0.07491) through the known relation 
(iV(123)|E? = i^frf ) |iV(123)) = 1(N\<t { *>t$\N), giving f nq = \U N . Because we neglect 
isospin breaking effects in this work, we chose the average pion mass m v = 138 MeV. 

C. Final states 

In Ref. || the wave function of the final NN-state was antisymmetrized and normalized on 
the nucleon level, assuming a plane wave with wave vector q in the relative coordinate r 
between the two nucleons. The coordinate representation of the nucleon-nucleon state vector 
I $aw(q)) was then written as 

(r|$Ar iV (q)) = $7V7v(q,r) = ^^=[ e ^ r -(-l)( s+T )e^ r ] ; 5 = 0,T = 1. (4) 

The full wave function of the final state respected the three-quark cluster nature of the 
nucleons 

(r | V NN (q, 123456)) = $^(q, r){iV(123)iV(456)} sr , (5) 

i.e. the nucleon wave function iV(123) was given by translationally- invariant shell-model 
configurations 



#(123) = \s 3 (b N ) [3] x , [1 3 ]c[21]t([21]ct)[21] 5 : [1 3 ]cts, LST= 0, 1/2, 1/2) T/5M 

= M123) • X C= ° ■ Xs7" ■ Xt7" ■ (6) 

S=— T= — 

X C= °i Xs z 2 arLC ^ Xt z 2 are color- singlet, spin and isospin three-quark states. $at(123) is 
the orbital part of the wave function, expressed in terms of the internal Jacobi coordinates 
p = ri - r 2 and A = r 3 - (ri + r 2 )/2 



$tv(123) = (v / 37r6^)~ 3/2 exp 



(7) 
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with a characteristic nucleoli oscillator parameter . This parameter does not have to be 
the same as the harmonic oscillator parameter b$ for the dibaryon wave function of Eq. ([!]), 
as it has been discussed in Refs. ||[7||. The relative Jacobi coordinate between the clusters 
is given by 

ri + r 2 + r 3 r 4 + r 5 + r 6 
r = . 8 

3 3 K J 

Note that the six-quark final state @ is antisymmetrized automatically when the state 
vector |^ArAr(q)) is substituted into the decay matrix element < ^AwCq); it \ O nq \ d'), because 
the initial state of Eq. ([I]) is fully antisymmetric. However, the antisymmetrizer-projector A 
(A 2 = A) contained in the initial state A \d!) = \d!) reduces considerably the normalization 
of the final state (it cuts all non-antisymmetrized parts of the cluster function of Eq. ([5]) 
which contain about 90% of the wave function - see below). Therefore, it is important 
to substitute from the beginning a final state wave function which is normalized (N) and 
antisymmetrized (^4) on the quark level. 

|^7VJv(q, 123456)) = MA{<& NN {q., r)iV(123)iV(456)}s T , ST = 01 . (9) 

The normalization factor M is determined by the standard orthonormalization condition 

(ttw(q')l*w(q)) = ( Wq')l Wq)> = <5 (3) (q' - q) , (10) 

which leads to 



({$ iViV (q / )iV(123)iV(456)}5T 


A\{$ NN (q)N(12S)N(456)}sT) 




$jviv(q)) 



The antisymmetrizer-projector A is 

^^(i-EE^rj^l 1 - 9 ^) ' A2=A - (12) 

P* CST is the pair-permutation operator for quarks i and j in orbital, color, spin and isospin 



space. It is instructive to calculate the normalization factor ( |Tl1) algebraically by factorization 



of the CST and X parts of the pair permutation P^ CST = P^ ST P<jf . The matrix element of 
P^ 6 ST between two NN-states in the ST=01 (or 10) channel is very well known (see e.g. flIB| ) 

({iV(123)iV(456)}5T=oi I Pg ST I {iV(123)iV(456)} 5T= oi) = ~ . (13) 

ol 

Inserting this value into Eq. (|TT[ ) reduces its right-hand-side to 

kt-2 = J_ 1 (^ J v(q / )^ J v(123)$ JV (456) | P£ | $ JVJ v(q)^(123)$ j y(456)) 

"!0 + 90 (M<rtlMq)> ' 

where $tv(123) is the orbital part of the nucleon wave function @ given in Eq. ([?]). The 
numerator in Eq. fll4D depends on the form of $jvjv(q, r), but for the plane wave in Eq. (f|) 
(or for any continuum wave function including f.s.i.) it is a finite value, i.e. it has to be zero 
compared with the ^-function in the denominator. Therefore, in our case the second term 
in Eq. fli"4]) vanishes and we obtain 

A/'=v / T0 (15) 

Note, that the expression (P£) = ( & NN $ N & N | P^ \ $ NN $ N $ N )/( & NN \ & NN ) receives its 
maximal value =1 in the special case of a Gaussian $Arjv(?") = (27r6^/3)~ 3//4 exp(— 3r 2 /46^). 
Therefore, for any relative NN wave function &nn we have the following constraints: 

< < 1 , or 9 < TV 2 < 10 (16) 



The value (|I5"D is equal to the usual identity factor J well-known in nuclear cluster physics 
(see e.g. [0]). It plays an important role for the projection of six-quark configurations onto 
baryon-baryon channels |20| . From now on we shall omit the antisymmetrizer A in front of 



the final state in the decay matrix element, but the identity factor ( |l~5|) may not be omitted: 

(d'\O wq (k)\V NN (<i),n) = (d'lO^VlOA^Mq)^) = VlO(d'\d nq (k)\^ NN (q),7r) (17) 

The inclusion of this factor, due to the antisymmetrization of the final two-nucleon wave 
function on the quark level, improves considerably the agreement of the results obtained in 
Ref. || with the experimentally suggested width. 
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D. Transition amplitude including intermediate states with up to two harmonic 

oscillator quanta 

As in Ref. |§, we calculate the decay matrix element of Eq. (|T7|) by inserting a complete 
set of six-quark configurations with quantum numbers of the final 1 S'o two-nucleon state 
(LST=001, J p =0+) 

(* NN (q),ir\6 irq (k)\d!) = VW y E <^ (q)Wl23)iV(456)} S2 y )1 |(n,6 6 ),{/},^T = 001) 

(«).{/} 

((n J 6 6 ),{/}, J LST=001|a w ,(k)|s^(6 6 )[51]x, [2 2 1 2 ] CT LST = 110, J p = 0") . (18) 

Here, {/} = {[/x]> L/cr]} an d (n) defines quark states with n harmonic oscillator (h.o.) 
excitation quanta, i.e. (n) = s 6 ~ n p n , s 6 ~ 2m (2s) m , (n = 2m), etc.; and is the h.o. parameter 
for the six-quark system. The summation in Eq. (|18D extends over a limited set of Young 
schemes [fx] and [far] '■ the possible representations of [fcr] m the sum in Eq. (|l8f) are given 
by the series of inner products of the [2 3 }c color and T=l [42] T isospin Young schemes 

[2 3 ] c o [42] T = [42] CT + [321] OT + [2 3 ] CT + [31 3 ] CT + [21 4 ] CT . (19) 

Only two spatial Young schemes [6}x and [42] x are compatible with the even-parity (L=0) 
N-N partial wave. Further constraints follow from the Pauli exclusion principle, i.e. [fx] ° 
[33] 5 o [fcr] = [l 6 ]xc5T- In the case of full spatial symmetry [6]x, only one color-isospin 
state [2 3 ]ct is allowed, but the Young scheme [42] x of the excited shell-model configurations 
is compatible with each state from the inner product given in Eq. (0). Our choice of a one- 
body transition (pion-production) operator defined in Eq. ([5]) further restricts the number 
of relevant intermediate states. The one-particle operator (|3|) can excite (or de-excite) only 
one quark of the initial s 5 p state. Therefore, the complete set of states in Eq. fll8f) is reduced 
to the configurations s 6 , s A p 2 and s 5 2s, knowing that higher one-particle excitations can be 
omitted because of a very small overlap with the final NN-state. Summarizing, the following 
intermediate states are taken into account in Eq. (|18|): 

• the energetically lowest (n=0) spatially symmetric state s 6 [6]x[2 3 ]cr, 

• the excited (n=2) translationally-invariant (orthogonalized to the 2S excitation of the 
six-quark cm.) state (s 4 p 2 — s 5 2s) with identical Young schemes [Q]x, [2 3 ]cr, an d 
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five excited (n=2) states s 4 p 2 [42]x[/cr] with CT Young schemes from the inner product 
of Eq. pi). 



It is interesting to note that all these configurations are also important for explaining 
the short-range nucleon-nucleon interaction. This was pointed out almost two decades ago 
[H,|r5,ZT| and thereafter discussed in many papers (see e.g. |2S| and references therein). 



Now we believe that a possible d! dibaryon has much potential for providing additional 
information on the innermost part of the nucleon-nucleon interaction, i.e. in the region 
where the nucleons overlap. 

E. Final state interaction 

To take into account the f.s.i. for the two outgoing nucleons, we consider separable-potential 
representations of the N-N interaction, namely the phenomenological potential of Tabakin 



22], and the separable model of Ueda et al. ||23|| , which is equivalent to the OBEP. The wave 
functions of the 1 Sq NN final states for the Tabakin potential are of the form 

®n~n(1i r) = (2tt) 3/2 cos(5 \ j (qr) - tan<5 n (gr) + A(q)- 



r 

-ar 



6 6 

+ BAq) cosar + B 2 (q) sinar} , (20) 

r r J 

while the separable potential model of Ueda et al. leads to the 1 5'o NN wave function 

-fr N p-PnT ' 



$^ (g,r) = (27r)- 3 / 2 cos5 \j ( qr )-t a n5oMqr) + A(q)^—3-J2Bn(q)^- \ ■ (21) 

n=l 



Here, 5o(q) is the phase shift of NN-scattering in the 1 So wave, and the functions A, A, 
Bi and Bi depend on the choice of parameters a, /?, 7 and $ for the two models (see the 
Appendix). 

We use the non-standard Tabakin potential because at short range, the NN wave functions 
obtained with this potential differ qualitatively from OBEP wave functions. In Fig. 0, the 
wave functions (E0) and (EI|) for both models are shown at an NN lab-energy of .5^^=100 
MeV. The relative wave function of Eq. (|20f) has a node at distances r w 0.4 — 0.5 fm (a stable 
position of the node in a large interval of NN-energies produces the same NN-scattering phase 



shifts as a repulsive core). The two models are phase equivalent, but differ in their off-shell 
behavior. In the following we will demonstrate that the results for differ considerably for 
both models, especially if the dibaryon mass comes close to the ttNN threshold. 



III. EXPLICIT CALCULATION USING THE FRACTIONAL PARENTAGE 
COEFFICIENT (F.P.C.) TECHNIQUE 



Our approximation for the decay amplitude in Eq. (|Tq) leads to a sum over products of two 
factors. The first factor is the so-called overlap integral of the intermediate six-quark config- 
uration with the outgoing two-nucleon state. The second factor is a shell-model transition 
matrix element that describes the production of the pion on a single quark in the dibaryon 
and the subsequent transition to an intermediate six-quark configuration. Both factors can 
be calculated with the standard fractional parentage coefficient (f.p.c.) technique, which was 
developed for quark-model calculations for example in Refs. |T3-p0f . 



A. Overlap integral of intermediate six-quark configurations with the NN-continuum 

In this subsection, we calculate the overlap integral of an intermediate six-quark configuration 
(n, bo) {/} with the (antisymmetrized and normalized) 1 S'o partial wave of the final NN-state 
introduced in Eq. (|j) 

<*^(?) I (n, &»){/})= (22) 
VlO($ L N = N (q)3{N((00,b N )l23)N((00, b N )456)} ST=01 | (n,b 6 ), [f x ], [f CT }, LST = 001) . 

Beginning with Eq. (^) we denote from now on the nucleon wave function of Eq. (^) of the 
translationally invariant shell-model as 

iV(123) = N({n'l' = 00, 6^)123) . (23) 

Here /' is the total orbital angular momentum contained in the internal Jacobi coordinates 
p, A, and b^ is the h.o. size parameter for the three-quark system. 

The overlap integral (^) is calculated using the standard f.p.c. technique for the quark 
shell model [ p1jpT| -|20| . For this purpose we use a f.p.c. decomposition of the six-quark 



in 



configuration (n, ^H/} into two three-quark clusters with CST quantum numbers of baryons 
{Bi((n'V, b 6 )123)B 2 ((n"l" , b e )456)}csT- Note, that for this procedure, the size parameter in 
the decomposition h§ differs from the nucleonic size parameter b^ 



\(n,b 6 )[f x ][f CT ]LST = 001)= £ £ J^[/«'^Kn") 



Bx(n')B 2 (n") 



x { ^~(r, fi/SI*) Y m {v) {B 1 ((n'l\ 6 6 )123) S a ((nT, 6 6 )456) W} (24) 



In expansion ([24]), </?j^r(r, //ro) is a h.o. wave function in the relative coordinate r of the 
two baryons with angular momentum L = L — (I' + Z") and iV = n — (n' + n") excitation 
quanta (/xro = ^2/3b 6 is the h.o. size parameter). As usual, rij x is the dimension of the 
IR [fx] °f the permutation symmetry group S$ for six particles p5] . rif> x and rij» are the 
dimensions of IRs [f' x ] and [/£] of the subgroups S' 3 and S3 in the reduction Sq D S3 x S3 . 
The coefficients C/g^ 2 and C}" } (n',n") are f.p.c. in the CST- and X-subspaces respectively. 
For simplicity, we omit in Eq. (|24D the indices for the dependence of U^^ 2 and C^(n',n") 
on the intermediate Young schemes /£, ST = f' x , f£ ST = f x , f' CT , fcTi f's, /t> /t> /c> 
and occuring for our chosen reduction chain of Eq. (0). 

With the help of (0), we can calculate the overlap (0). The three-quark - three-quark 
decomposition is of course the most adequate expansion for projecting onto the NN-channel. 
The projection for a given intermediate state (n){/} 

$f n =y { ° /} (r) = v / 10({iV((00,^)123)iV((00,6^)456)} 5T=01 1 (n,b 6 )[f x ][f CT ]LST = 001) (25) 

receives non-vanishing contributions only from NN-components in Eq. (|24| ) because non- 
nucleonic clusterings, such as -E>i(123)-E>2(456), are orthogonal to iV(123)iV(456) in CST 
space. Furthermore, the overlap of excited nucleonic clusters, e.g. iV((20, 6 6 )123), with the 
ground state nucleon iV((00, 6^)123) can be neglected if we assume that the size parameter 
&6 of the six-quark configuration (|4]) does not differ considerably from the quark core radius 
b jv of the nucleon. In fact, because be 7^ bjy, the non-zero overlap integral between excited 
and non-excited nucleons is 



(bi/b 2 N -l) ( 2b 6 /b N 
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(Af((20, f> 6 )123) I Af((00, f>„)123)) = jf^JJ^j ( J ' < 26) 



1 I 



The sum over all possible terms gives a negligible contribution to the final result because 
the different terms interfere destructively (see next section). Due to these restrictions we are 
lead to the expression 



$(nX/}( r ) ~ (^((00,6^)123) | iV((00A)123)> (N((00, b N )A5Q) \ N((00, 6 6 )456)) 

x y/wJ-Lu$ Cg(0, 0) Vn0 (r, ^2/3b 6 ) Y 00 (r) (27) 

V n fx 



where 



(0,0) 



1, if n = 0, [f x ] = [6] 
% n = 2, [f x ] = [6] 



n 



2, [fx] = [42] 



The coefficients C^(n',n") are calculated by general methods from the TISM (see, e.g. Ref. 



. The values of are given in Table [I]. The general rule for calculating f.p.c. in the 
CST-subspace is the factorization of the value Uy^ 2 fT4}{i~8f (symbolically) 

Ufff 2 = SFct»s SFc,t (28) 

in terms of scalar factors SFct»s and SFq,t of Clebsch-Gordan coefficients of the unitary 
groups SU(12)cst and SU(6) C t for the reductions SU(12)cst 3 SU(<o) CT x SU(2) S and 
SU(6) ct 3 SU(3)c x SU(2)t, respectively (which are links of the common reduction chain 
(0)). The necessary SFs are tabulated in Refs. p5|-p0|. With expression (0), the overlap 
integral of Eq. (|22|) reduces to 



<*r»°(9) i («,*.){/}) = <*r»°(9) i *(s«> 

^(ifi^) 8 (^f^M • (20) 

Eq. (|29|) contains a simple radial integral 

/8^\ 3/4 



4t (q) = V^ r 2 dr $^°(g, r> n0 (r, y^fee) , (30) 



V 3 / 

which can be calculated analytically for a plane wave $^*(g,r) = (27r)~ 3 / 2 jo (<??"), as well as 
for the Tabakin f.s.i. wave functions of Eq. ( |20"D and for the Ueda f.s.i. wave functions given 
in Eq. (^1|). Results for Initio) are listed in the appendix. The large bracket in Eq. (|29|) 
involving the ratio of the two h.o. size parameters b^/bg comes from the overlap of the two 
nucleon clusters (iV((00, 6jv)123) | iV((00, 6 6 )123)) (JV((00, b N )A56) \ N((00, 6 6 )456)). 

1? 



B. Shell- model transition matrix element 



The shell-model matrix element of the pion-production operator O nq defined in Eq. (0) (the 
second factor in the decay amplitude introduced in Eq. (0)) is proportional to the one- 
particle matrix element of the spin-isospin-flip operator a^Tj K \ The remaining five quarks 
act as spectators for the transition. Chosing the sixth quark, we write the matrix element 
for the emission of a 7r~: 

((n, b 6 ){f}LSf = 001, J p = 0+ T z = l; 7r~ | a^(k) | d') 

= 6— i k ((n,b 6 ){f}LSf = 001,T z = l|( O - 6 -k)r 6 (+) exp^k-p 6 ) | d!) (31) 
m n v /2^ 7r (27r) 3 V 6 / 

Here, p 6 = r 6 — | Y%=i v ii is the Jacobi coordinate, and <x 6 and r 6 are spin and isospin of 
the sixth quark. The momentum of the pion is k = fck, and the factor 6 in front of the 
one-particle matrix element contains the summation over all six quarks. The natural choice 
for a f.p.c. decomposition is clearly the separation of the last quark g 6 — > q 5 x q (one-particle 
f.p.c), which allows to exploit the orthogonality constraints for the five spectator quarks. 
With the one-particle f.p.c. expansion of the shell-model states, the right-hand side of Eq. 
(|31|) reduces to a sum of one-particle spin-isospin-flip amplitudes with algebraic coefficients: 



({n,b 6 ){f}LSf = 001, T z = l | (<r 6 ■ k)ri +) e- l i^ \ d') = £ £ ^ 



f' x V U fx f' x V U fx 

x £ U {f - } (S'f', S 6 = T 6 = 1/2) £ U {f} (S'T', S 6 = T 6 = 1/2) £(1, —M, 1,M\0, 0) 

S',f' S",T' M 

x {5„, 2 4^(^,5) Cg ) ( S 4 p,p)X 6 (k;llM,00) 

+C^ 1]x (s 5 ,p) [S n , C$ x (s 5 , s) + 5 n , 2 C$ x (s 5 , 2s)] X 6 (k; 11M, n0)} 
x3i(S',S e : 5 = | (cr 6 ■ k) | 5", 5*6 : S=l,M){f',T 6 :f = T z = l\ r 6 (+) | T\ T 6 : T = 0) (32) 

where the functions Xg(k; 11M, nO) are spatial integrals 

X 6 (k;llM,nO)= / (k ■ p fl )y 1 _ M (p 6 ) d 2 p 6 

x / pldp 6 (p n0 (p 6 , y / 6/56 6 )v9ii(p 6 , \fo/5b 6 ) ji(5/6fcp 6 ) (33) 



Here 5" and T" (5*' and T') are spin and isospin of the five spectator quarks after the 
separation of the sixth quark. The first term in curly brackets on the right- hand- side of 
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Eq. ( p2[) corresponds to the quark transition from an initial s state to a final p state. The 
second term corresponds to the quark transition from an initial p state to a final s or 2s 
state. In Eq. ([52]) , we use almost the same notation for one-particle f.p.c. U^(S'T', SqTq) and 
Cfe((n'), (n")) as before for the three-particle f.p.c. 's in Eq. (|24|). Note, that the one-particle 
f.p.c. L r {/}(S"T', SqTq) in the CST-subspace can be calculated with one-particle scalar factors 
SF (cf. Eq. (|28|) ), which can be found for example in Ref. [19| . 



Due to the orthogonality restrictions for the five spectator quarks, the summations over 
fx,f' x ,S',T',S',T' collapse to Sf j> Ss',s' ^s',s 6 $f',T' ^f',T 6 y an d the only non-vanishing 
elementary spin- and isospin-flip amplitudes in Eq. (|32"D are 

(S' = S 6 = 1/2 : S = | \S' = S 6 = 1/2:S=1, —M) = -(-1)" 5^ 
(f = T 6 = 1/2 : f = 1,T Z | r 6 W | T = T 6 = 1/2 : T = 0) = (-1) K 6 K * (34) 



C. Decay amplitude after summation over allowed intermediate states 

Collecting the shell- model matrix element of Eqs. (|3T| ) -(|32| ) for pion production and the 
overlap integrals of Eq. (|29|) for all intermediate six-quark configurations, and performing 
the remaining radial integrals in Eqs. ( p0|) and (|33D, leads to the following result for the full 
decay amplitude defined in Eq. (|T8|) 



<*jvjv(q),7r-|a 9 (k)|d') 



lOi 



27(2tt) 3 / 2 V3tt 



3/4 




f-rrq 

m T V E,, 



266/6. 



iV 



1 + 6 2 



(A;6 6 ) 2 exp(-^^ 



6 / N , 



24 



x 



AO) 

NN 



(?) + 



27 V 24 

The overlap integrals Inn(q) an d Inn(q) can be found in the appendix. Note that the 
inclusion of overlap terms with excited nucleon configurations in the intermediate six-quark 
states, originating from the different harmonic oscillator parameters b§ ^ 6jy, given in Eq. 
(|26| ) leads to a non-essential renormalization factor of the decay amplitude (|3 



(35) 



576(1 + 61/6?,) V 6 fc ° 6 



in front of the term I^ N (q) on the right hand side of Eq. 



(36) 



This factor (|36| ) can be 



omitted for small kbg. The k 2 behavior of the decay amplitude is due to (i) a factor k in the 
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transition opertor of eq.(|3]), (ii) due to the fact that a p-wave quark is involved either in the 
initial or the final state of the one-particle transition matrix element on the right-hand-side 
of the Eq. fl32|). We recall that in the case of the A-isobar decay into the irN channel, the 
transition matrix element is proportional only to k 1 , corresponding to the <Tj ■ k term in Eq. 
(|3|). The k 2 behavior of the d! decay amplitude (^) leads to a very strong dependence of 
the d! decay width on the value of Mai, as we will see in the next section. 



IV. NUMERICAL RESULTS AND DISCUSSION 



The total hadronic decay width of the possible d! dibaryon IV contains three partial widths 



^nn ST,,-— pp 



(37) 



which are equal to each other IV 



/</< 



^ir°pn 



when we neglect isospin breaking 



effects. The partial it pp decay width IV-™, is defined by the standard expression 



IV PP = 2tt J d 3 q J d 3 k5lM d , - 2M N 



k 2 



AM 



N 



M, 



N 



m 2 + k 2 



x 



(*jvjv(q),7r-|a 9 (k)|d') 



(3f 



where q = qi 2 q2 is the relative momentum of the two final protons and k is the momentum 
of emitted pion in the cm. of the d! dibaryon. The 5-function conserves the energy in the 
decay, while the integration over the momentum conserving 5^(qi + q2 + k) has already 
been exploited [0] in Eq. The integration over three-particle phase space leads to the 

following result for the partial d' — > ir~pp decay width 



2 5 10 2 v / 6A 1 ( 2b 6 /b 
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A' 



7T pp 



3 8 ^ 4n ml \l + b 2 6 /b 



N , 



'I ma x 



2M N (k b e 



exp 



„ 2M N + Jml + k 2 * V 12 



k 2 h 2 



q 2 dq. (39) 



Here, energy conservation relates the pion momentum ko to the NN relative momentum q 



via 



AM, 



N 



Mai 



M, 



N 



\ 



Mai 



M 



Mai - 2M, 



N 



N , 



M 



+ mi 



N , 



1/2 



and for q max = ^M^{Md' — 2Mjy — rn n ) } all available decay energy is converted to kinetic 
energy in the relative NN-system, and none to the pion E n = m n , ko = 0. 

The calculated decay widths are shown in Table O, where we have introduced the abbre- 
viations p.w., T and U. Here, p.w. refers to a calculation employing a plane- wave final N-N 
state (Hj), while T and U refer to calculations using the Tabakin ||] (T) and Ueda et al. g| 
(U) separable NN potentials for the final state interaction. 

In parentheses we give the results obtained in the approximation of using only one in- 
termediate six-quark configuration s 6 (n=0). With the exception of the results for the Ueda 
NN-potential for sets 1, 2 and 4 (for which the d' mass is 400 - 650 MeV above the 7rNN 
threshold), the inclusion of all Pauli principle allowed intermediate 2hu shell-model configu- 
rations tends to increase the decay width by some 20 - 30 %. The largest effect is obtained 
for d! masses rather close to threshold, exemplarily shown for sets 3 and 5. 

It can be seen from Table [TI] and Fig. [| that the pionic decay width of the d' is very 
sensitive to the dibaryon mass M^, which determines the available phase space of the three- 
body tcNN decay. The sensitivity grows dramatically near the irNN threshold (2016 MeV). 
If we extrapolate the results of Table |T| to the experimental value of M# = 2065 MeV, we 
obtain a very strong reduction of IV as compared with the quite realistic variants (sets 3 
and 5) in Table yj: 

T p d r = 0.032 MeV, T T d , = 0.046 MeV, Y u d , = 0.083 MeV, if 6^=0.595 fm and 6 6 =0.95 fm 
T p d , w - = 0.018 MeV, T% = 0.045 MeV, 1^ = 0.040 MeV, if 6^=0.6 fm and 6 6 =1.24 fm . 
This strong dependence of on the value of is a consequence of the high power of 
(&; &6) 5 m the integrand of Eq. d39|) . The origin of this k\ behavior (compared with a k\ 
behavior in case of the A-isobar decay) comes, as explained above, from the necessity to 
excite (or de-excite) a p-wave quark for the production of a pion. Note that for small q ma x 



(when is close to the 7rNN threshold) the function ko(q) is linear in the factor Jq^ax ~~ Q 2 



and can be written as k (q) m q m ax\J 4771^(1 — q 2 / q^^) / M d i . Therefore, for small q max the 
integral in Eq. ( |59|) behaves as q^ax- The second high-power factor in Eq. (|39|) is the scale 
factor 

/ 2h/b N A 1 ' 



l + b 2 



6 / N , 



which depends sensitively on the ratio b%/bj^. However, this sensitivity is considerably re- 
duced by the factor 6g in the integrand. The product 

5 / 2b 6 /b N y 2 
6 {i+bi/b%) 

is a quite smooth function of b^/b 6 . For 6^=0.6 fm this product varies from 0.078 fm 5 to 
0.158 fm 5 , if &6 varies from 0.6 fm to 1.24 fm. 

For small q max , f.s.i. make an important contribution to the d! decay width because of 
the large scattering length in the 1 5'o wave a s = —23.7 fm. The f.s.i. enhances the decay 
width for example by about 85% for set 5 in Table [TI|. At the experimental mass M^=2065 
MeV, the hadronic decay width is more than doubled by the final state interaction. It is 
interesting that in the case of the Tabakin model with a nodal NN wave function at short 
range, the contribution from f.s.i. is smaller than for the Ueda model and can even decrease 
the width compared to the plane wave result (cf. set 4). This is a direct consequence 
of an approximate orthogonality of the nodal wave function of the Tabakin model to the 
projection of the intermediate s 6 configuration (i.e. the h.o. function y? o) of Eq. (|27|) onto 
the NN channel. This can easily be seen from Fig. |2|, where both wave functions are shown. 
The approximate orthogonality of the functions ipoo and $™ afcm in the integrand of Eq. 
(29) reduces considerably the overlap factor I^ N (q), which gives the dominant contribution 
to the d! decay width (see values in parenthesis in Table 0). As it can be seen in Fig. |3], 
the disagreement between the Tabakin and Ueda models grows with increasing dibaryon 
mass Ma' (the influence of the large scattering length a s , which is common for both models, 
becomes negligible compared to the effect of the larger phase space). For sets 1, 2 and 4 in 
Table ||, the Tabakin model leads again to values of IV, which are even smaller than IV in 
the plane wave approximation neglecting f.s.i. 



V. SUMMARY 

In this work we have studied the pionic decay of a possible d' dibaryon within the micro- 
scopic quark shell-model. We use a single-quark transition operator which describes the 
production of the pion on a single quark. The dibaryon wave function is given as a single 
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six-quark translationally invariant shell model configuration, which has been found to pro- 
vide an adequate description of the d' Previous results from H have been improved 
mainly in three points, leading to a complete calculation in the sense, that i) the calculation 
is performed consistently on the quark level, i.e. the final two-nucleon state is normalized 
and antisymmetrized on the quark level, ii) all important intermediate six-quark states with 
nonvanishing overlaps with the final two nucleons are included, and iii) the strong final state 
interaction for the two nucleons is taken into account on the basis of the separable Tabakin 



22| and Ueda et al. |23[ NN potentials. 

Not surprisingly, the small available phase space in the three-body decay is the dominant 
mechanism for the narrow width of the d' . The large identity factor (15), on the other 
hand, enhances the results of previous evaluations disregarding the identity of quarks from 
different nucleons. The inclusion of all Pauli principle allowed intermediate 2hu shell-model 
configurations tends to increase the decay width by some 20 - 30 %. Furthermore, final 
state interaction for the two outgoing nucleons also increases the decay width considerably, 
if the d! mass is close to the ttNN threshold. Due to these three effects, the calculated pionic 
decay widths lie between = 0.18-0.32 MeV for the most realistic set 5, having a d' mass 
close to the experimentally suggested one and a characteristic hadronic size of the dibaryon 
of &6 ~ 1 f m - This qualitatively agrees with the experimentally suggested value Yd> = 0.5 
MeV. 

Despite the fact that both, the Tabakin and Ueda f.s.i. models, are unsuitable for large NN 
energies (as in parameter sets 1, 2 and 4), the two models demonstrate the strong influence of 
the short-range behavior of the NN wave function on the d! decay width. (See e. g. Figs. [|and 
U). Recall that these two models are typical representatives of qualitatively different classes 
of NN phenomenology. Whereas the Ueda separable potential is an approximation of the 
OBEP, i. e. a model with short-range repulsion, the Tabakin potential can be considered as 
a unitary-pole approximation (UPA) [|2fJ of a Moscow-type potential model |27j with short- 
range attraction and forbidden states. The Moscow model proceeds from the assumption 
of a six-quark origin of the short-range NN interaction and pretends to give an adequate 
description of the non-local character of the NN force. The main conclusion to be drawn 
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here is, that these two models, which are phase equivalent, differ considerably in their effect 
on the d! decay width. Therefore, a possible d! dibaryon would provide a natural laboratory 
for detailed studies of the short-range NN interaction. 

An interesting continuation of this work would be to go beyond phenomenological NN 
potential models and use a completely microscopic quark model approach (see e.g. [f2^] and 
references therein). For example, one could calculate the pionic decay of the d' dibaryon 
using a final NN-scattering wave function that is based on the same microscopic quark 
Hamiltonian which simultaneously describes the mass M<# and structure of the d' dibaryon. 
However, such a calculation is complicated by the fact that we have used two different 
Hamiltonians, i. e. two different confinement strengths, for three-quark baryons and the six- 
quark d! dibaryon . Thus, the d! dibaryon could not be explained in terms of the standard 
constituent quark model, using a common Hamiltonian for any number of quarks. On the 
other hand, if the d! really exists, this may be taken as an indication that the effective 
(nonperturbative) quark-quark interaction depends on the state of the system. 

APPENDIX A: 

In this appendix we present the analytical expressions for the radial integrals I^n^) defined 
in Eq. (|3(]), which are needed to calculate the overlap integral of Eq. (|22|) between differ- 
ent intermediate six-quark shell-model configurations and the two outgoing nucleons. We 
recall, that the relative NN-wave function may be described by a simple plane wave (p.w.) 

r ) = (27r)~ 3 ^ 2 io(?^) or by f.s.i. wave functions resulting, for example, from separa- 
ble potential representations of the NN-interaction. We use the results for the NN-wave 
function obtained by Tabakin [22] given in Eq. (p0|) and the result obtained by Ueda and 



co-workers ]23[ given in Eq. (^1|). Note that coefficients A and B n in Eq. ( pT|) depend on 



parameters of the one-term separable potential of Ueda et al. |23| V(q, q') = —Mug(q)g(q'), 

(<? 2 +/3 2 )(<? 2 +7 2 )' 



( 2 2\ ~ ~ 

9((l) = J2 n ( 2?a2u~2i 2\ - F° r A and B n we use the following expressions: 



a(\ = 2 9(g) \- Cn ( l 2 + qV 
{q> a 1-G(g)^7 2 -^ W + q\ 



1Q 



B n (q) = a 



2 g(q) 



'PI + q, 



where a 2 = 2-K 2 m N M 11 /% 2 , G{q) = \a 2 f™[{g 2 {k)k 2 - g 2 {q)q 2 )/(k 2 - q 2 )]dk and a value of 
7, which is not fixed in the ref. p3j , is fitted to the singlet scattering length (a s =-23.7 fm), 
7 = 11.114 fm" 1 . 



In Table III we introduced the following abbreviations: 



gV>(x) 



(0) 



(2) 



-x 2 /3 



^ - e -?l*Im*{£ 



X\/TC 



VV3, 



\fl (I 4 



x 



XWTC 



X^/ft \3 9 

e - 2 /3 [i_$( x / 3 )] 



VS (\ 4 



2 \xyfH V3 9 



- + -x' 



1 _ !^ e ^ 2 /3 /m$ (™\ 



! + £/) e ,V3 [i_$( x /3)]|, 



+ (sin x 2 /3 - cos x 2 /3) [l - i?e$ ((1 + z)x/V3 



2 (x\/7r 

x 2 ^) sin x 2 /3 + ( — I — x 2 ) cos x 2 /3 

3 9/ 1 V3 9 1 

/^/- 4 x 2 ^- 
'2| 9 2h/tt 



(sin x 2 /3 + cos x 2 /3j Jm$ ((1 + %)xj , 
(sin x 2 /3 + cos x 2 /3) [l - i?e$ ((1 + i)x/Vs)] + 
(sin x 2 /3 — cos x 2 /3^ Jm$ ((1 + i)x/ a/3) , 

° ^ X ^ ^ + j^x 2 ) sin x 2 /3 - f J - j^x 2 ) cos x 2 /3l [l - i?e$ ((1 + i)x/V3 
39/ V3 9 / - 

Jm$ ((1 + i)x/V3 

- — -x 2 ) sin x 2 /3 + I - — -x 2 ) cos x 2 /3 1 - i?e$ ((l + z>/V3 
3 9 / V3 9 / J l ^ 

^ + ^x 2 ^j sin x 2 /3 - Q - ^x 2 ^) cos x 2 /3 im$ ((1 + i)x/v / 3) | 
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I. T. Obukhovsky et ai., Figure 1, 

"Pionic decay of a possible d! dibaryon and ..." 




d' 

FIG. 1. Quark line diagram of the pionic dibaryon decay. The elementary pion is produced on 
a single quark, leaving the remaining six quarks in a relative 1 Sq nucleon-nucleon scattering state. 



I. T. Obukhovsky et ai., Figure 2, 

"Pionic decay of a possible d! dibaryon and ..." 
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FIG. 2. Wave functions of the final 1 Sq state for two NN interaction models [22,23] at fixed 
lab-energy £tv7V=100 MeV. The projection of the s 6 six-quark configuration onto the NN channel 
is also shown. 
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I. T. Obukhovsky et ai., Figure 3, 

"Pionic decay of a possible d! dibaryon and ..." 
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FIG. 3. Pionic decay width IV of the d! as a function of the dibaryon mass for different final 
state interactions (f.s.i.) between the outgoing nucleons: (i) plane wave (p.w., dotted curve), (ii) 
with f.s.i. using the Tabakin [22] potential (dashed curve), (iii) with f.s.i. using the Ueda et al. [23] 
potential (plain curve). The harmonic oscillator parameters 6^=0.595 fm and &6=0.95 fm are those 
of set 5 in Table 2. The 7rNN threshold for the decay is at 2016 MeV, while the experimentally 
suggested resonance position of the d! is at 2065 MeV. 
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Tables: 



[fcrs] 


[1%TS 


[2 2 1 2 ]cts 


[for] 


[2 3 ]ct 


[42] ct 


[321] CT 


[2 3 ]ct 


[31 3 ]ct 


[21 4 ]ct 


ttNN 
U {f} 


fl 
V 9 


/IT 

V 20 


/l6 
V 45 




FT 

V 18 






TABLE I. The CST part of the f.p.c. three-quark - three-quark decomposition for the projection 
onto the NN-channel UM , {/} = {[fcrs], [fcr]} 







IV pp [MeV] 


fV[MeV] 


set 


b N [fm] 


hifm] 


M d /[MeV] 


p.w. 


T 


U 


p.w. 


T 


U 


1 


0.45 


0.59 


2705 


56.8(44.1) 


7.8(5.2) 


41.2(46.0) 


170.5 


23.3 


123.6 


2 


0.47 


0.65 


2680 


44.2(35.7) 


8.3(7.0) 


32.6(36.5) 


132.5 


24.9 


97.8 


3 


0.6 


1.24 


2162 


0.22(0.17) 


0.27(0.18) 


0.25(0.22) 


0.67 


0.81 


0.76 


4 


0.595 


0.78 


2484 


28.3(22.6) 


10.9(8.9) 


21.8(21.9) 


84.8 


32.6 


65.4 


5 


0.595 


0.95 


2092 


0.058(0.036) 


0.061(0.049) 


0.107(0.071) 


0.173 


0.183 


0.321 



TABLE II. Calculated 7r~ decay width T^- pp and the total hadronic decay width of the d' 
dibaryon for five different d! masses and wave functions (b§ is the characteristic d' size parameter, 
and is the quark core radius of the nucleon). Masses and wave functions of the d' were obtained 
in Refs. [6,7] within different models for the microscopic q-q interaction. 



9.7 



Model 


4Ui) = 


4 2 k?) = 


p.w. 




f (2) (qbe) 


Tabakin 
Eq. (0) 


cos5 [f^(qb 6 )+tanS g^(qb 6 )+ 
A(q)aF^(ab 6 ) + B x {q)pG { ? ) (/36 6 ) + 
B 2 (q){3G { 2 °\(3b 6 ) 


cos5 [f ( - 2 \qb 6 ) + tan5 g^\qb (i )+ 
A(q)aF^(ab 6 )+B 1 (q)PG^ ) (Pb 6 ) + 
B 2 (q)(3G^(/3b 6 ) 


Ueda 
Eq. (H) 


cos5 [f^(qb 6 )+tanS g(°\qb 6 )+ 
l(g) 7 F(°)( 7 6 6 ) - E n Bn(q)(3nF^((3 n b 6 ) 


cos5 [f^(qbe) + tan5 g {2 \qbe)+ 
Aiq^FW^h) -J2 n B n (q)f3 n FW(p n b 6 ) 



TABLE III. Radial integrals I^(q) for pi ane waves (p.w.) and f.s.i. functions given in Eqs. 



(20) and (21) for the Tabakin [22] and Ueda [23] separable potential formulation of the relative 
NN-wave function. 



